%fractionS0 is the vector ratio between S_t and S_0
%time_interval is the time we want to interpolate
%%
function sigma = estimateSigma(fractionS0,time,volatility,St,time_interval)
    %Interpolate sigma
    sigma = interp2Optimized(fractionS0,time,volatility,St,time_interval,'spline');